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We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid 
spacetimes with a helical Killing vector k" , heuristically, the generator of time-translations in a 
corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have 
asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black- 
hole binaries, a rigidity theorem implies that the Killing vector lies along the horizon's generators, 
and from this one can deduce the zeroth law (constant surface gravity of the horizon). Remarkably, 
although the mass and angular momentum of such a system are not defined, there is an exact first 
law, relating the change in the asymptotic Noether charge to the changes in the vorticity, baryon 
mass, and entropy of the fiuid, and in the area of black holes. 

Binary systems with AIQ small have an approximate asymptopia in which one can write the first 
law in terms of the asymptotic mass and angular momentum. Asymptotic fiatness is precise in 
two classes of solutions used to model binary systems: spacetimes satisfying the post-Newtonian 
equations, and solutions to a modified set of field equations that have a spatially conformally flat 
metric. (The spatial conformal flatness formalism with helical symmetry, however, is consistent 
with maximal slicing only if replaces the extrinsic curvature in the field equations by an artificially 
tracefree expression in terms of the shift vector.) For these spacetimes, nearby equilibria whose stars 
have the same vorticity obey the relation 5M = Q5J, from which one can obtain a turning point 
criterion that governs the stability of orbits. 



I. INTRODUCTION 

Beginning with papers by Blackburn and Detweiler [|l|j2| , several authors have used spacetimes with a helical Killing 
vector 1^ to model binary systems in the context of general relativity. Such spacetimes can be regarded as having equal 
amounts of incoming and outgoing radiation; they are a counterpart in general relativity of the stationary solution 
due to Schild Q| that describes two oppositely charged particles whose electromagnetic field is the half-advanced -I- 
half-retarded solution of the orbiting charges. Because the radiation field of such a stationary solution has infinite 
energy, spacetimes that describe the corresponding general relativistic binaries are not asymptotically flat. Instead, 
the asymptotic mass rises linearly with a naturally defined radial coordinate. 

The formal lack of asymptotic flatness has been handled in several related ways. As Detweiler has emphasized, 
one can deflne an approximate asymptotic region for systems in which the energy emitted in gravitational waves in 
a dynamical time is small compared to the mass of the system. In this "local wave zone," the geometry describes 
gravitational waves propagating on a Schwarzschild background. In the more restrictive context of the post-Newtonian 
approximation, one regains asymptotic flatness, because there is no radiation through second post-Newtonian order. 
Finally, a number of authors considered spacetimes with conformally flat spacelike slices that satisfy a truncated 
set of field equations consisting of the initial value equations, a field equation for the spatial conformal factor, and 
the equation of hydrostatic equilibrium. Like the post-Newtonian spacetimes, these conformally flat spacetimes are 
nonradiative and (as we show) are asymptotically flat. 

We consider binary systems modeled in the exact theory (without asymptotic flatness) and then apply our results 
to the post-Newtonian spacetimes and spatially conformally flat spacetimes that retain asymptotic flatness. In each 
case one uses the Killing vector, /c", to define a conserved current and associated charges. For the exact vacuum and 
perfect-fluid spacetimes, the Noether current of the helical Killing vector assigns to each spacetime a charge Q. (See, 



^In a spacetime with a rotational Killing vector 0" and a timelike Killing vector t", each combination t" -\- flcj)" , with fl 
constant and nonzero, will be called a helical (or helicoidal) Killing vector. We give a precise definition in ^ and discuss its 
relation to a previous definition by Bonazzola et al M 
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for example, Refs. |10-15|.) Despite the lack of asymptotic flatness one can choose the current to make Q finite, and 
it Q is independent of the 2-surface S on which it is evaluated, as long as S lies outside the matter and all black 
holes. The Noether current assigns to each black hole a charge that can be identified with its entropy (its area, in the 
spacetimes we consider); and we obtain a version of the first law (Eq. ( |55| ) below) that expresses the change 5Q in 
terms of changes in the vorticit^ baryon mass, and entropy of the fluid, and in the area of black holes. Independent 
work by Baker and Detweiler |16|] obtains a similar first law for spacetimes with approximate asymptotic flatness at 
finite distance from the binary. 

In the asymptotically flat spacetimes mentioned above, the helical Killing vector has the asymptotic form 
/c" = + where t" and (jf generate asymptotic symmetries associated with time-translations and rotations. 

Neighboring perfect-fluid equilibria in a post-Newtonian or a spatially conformally flat framework satisfy a first law 
of thermodynamics of the form 

5M = nSJ + I \fMS + flAdMB + v^AdCJ + V —k^SA,. 

Here M and J are the ADM mass and angular momentum of the spacetime (see Eqs. ( |108| , |109D ; T and p, are the 
redshifted temperature and chemical potential; (IMb is the baryon mass of a fluid element; and dCa is related to the 
circulation of a fluid element (see Eqs. ( [5^j57| )). 

Note that, in the full theory, models of binaries with a helical Killing vector can only have corotating black holes. If 
their generators do not lie along the Killing vector the black holes will have nonzero shear and thus (assuming positive 
energy) increasing area; and this is inconsistent with the assumption of a helical Killing vector. In an appendix, we 
derive a virial relation for binary neutron-star systems in a conformally flat framework and show that the relation is 
equivalent to the equality of the Komar and ADM mass. 

One other class of asymptotically flat spacetimes with a single Killing vector is worth mentioning. These are non- 
axisymmetric stars whose figure is stationary in an inertial frame, the analog in general relativity of the Newtonian 
Dedekind ellipsoids. We expect that such stationary, nonaxisymmetric perfect-fiuid spacetimes exist; their velocity 
fields have nonzero shear, however, and cannot be stationary when viscosity is present. ]l7[ | 

Conventions: Spacetime indices are Greek, spatial indices Latin, and the metric signature is h ++. Readers 

familiar with abstract indices can regard indices early in the alphabet as abstract, while i,j, k, I are concrete, associated 
with a chart {a;*}. We use the dual form of Stokes' theorem for the divergence of an antisymmetric tensor A'^'"^'^ , 
namely 



[ V^A°'-^''dSo....f3 = / A" 



■f^^'dS, 



Q---/37 



where dSa = ^aPjsdS^'*^ , dSap = Cafs-ysdS^^ . For example, in an oriented chart t, {a;*} with S a surface of constant t 
and dT, a surface of constant t and x^, dSa — at\/—gd^x, dSap — -^iS atV px^ — V aX^V pt)^ — gd^ x . Finally, if S is 
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a 2-surface in a 3-space S and eabc is the volume form on S associated with a 3-metric 7afc, we write dSa = eabcdS^' 
for S a surface of constant r, dSa = "S/ ar^/^d^x. 



II. HELICAL KILLING VECTORS, EVENT HORIZONS, AND THE ZEROTH LAW 

We consider globally hyperbolic spacetimes Af, gap that have a symmetry vector fc", a Killing vector that generates 
a symmetry of the matter fields. Our particular interest is in stationary binary systems, systems whose Killing vector 
fc" has helical integral curves with a fixed period T; but our results hold for a broader class of spacetimes with a single 
Killing vector. 

We begin by using the periodic orbits just mentioned to define a helical vector field. We want a definition that agrees, 
for stationary, axisymmetric spacetimes, with Killing vectors of the form t" + fit/)", where t"' is the asymptotically 
timelike Killing vector and 0" the rotational Killing vector. Let Xt be the family of diffeos generated by fc", moving 
each point P g M a parameter distance t along the integral curve of fc" through P. Although a helical vector is 
spacelike at distances from the axis larger than T/27r, its integral curves spiral each period to points that are timelike 
separated from their starting points; at least they are timelike separated when one is outside a finite region that 
encloses any horizon or ergosphere. Without this last caveat, one could define a helical vector by the requirement 
that, for each point P, Xt{P) be timelike separated from P. To include the caveat, one requires that the condition 
hold only outside some sphere. Let 5 be a spacelike sphere, and let T be the timelike surface swept out by the action 
of Xt on S: T{S) = UtXt{S); we call T the history of S. 
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Definition II. 1 A vector field k" is helical if there is a smallest T > for which P and Xt{P) ore timelike separated 
for every P outside the history T of some sphere. 

When the spacetime admits a foliation by timchke hnes, this definition is equivalent to the following definition, 
essentially that of Bonazzola et al. 

Proposition II. 1 A vector field k" is helical if it can be written in the form 

k°' = t^" + , (f) 

where is spacelike and has circular orbits with parameter length 27r, except where it vanishes; n is a constant; and, 
t" is timelike outside the history T of some sphere. Conversely, if a vector field is helical, and if the .spacetime can be 
foliated by timelike curves that respect the action Xt{P), then k" can be written in the form ^ 

Because there are spacetimes with helical vectors that do not allow foliations respecting the action xt{P), the 
Bonazzola et al. definition is slightly more restrictive than ours; they are also more restrictive in requiring the exis- 
tence of a 2-dimensional submanifold, the axis of symmetry, on which 0" vanishes; and in requiring that t" be timelike 
everywhere. Note that, although the proposition displays the intuitive character of a helical vector, t" and are far 
from unique. Each foliation of M by a family of timelike curves that respects xt gives a different decomposition of 
k" of the form (§). 

Proof of Proposition. The first part of the Proposition, that a vector of the form is helical, is immediate. We 
prove as follows that a helical vector can be written in this form. Define a scalar by requiring it to have the value t on 
Xt(S]), with S a Cauchy surface. Let be the vector tangent to our congruence of timelike curves, each parametrized 
by t. 

Let ipt be the family of diffeos generated by P. Each integral curve of fc" can be projected to a circle on S by 
pushing it down to S along the timelike congruence: The circle through each point P G S is 

t^c{t) ■.= ^_toxt{P)- 

One obtains a circle with parameter length 27r by reparametrizing c, defining C(s) :— c[sT/(27r)]. Finally, define 0" 
on S as the vector field tangent at each point P to the circles C{t) through P; and drag 0" by V't to extend it to M. 
Then k°' ^ t" + Qcj)'^, with ft = 2n/T. □ 

In particular a spacetime that is stationary and axisymmetric, with asymptotically timelike Killing vector t" and 
rotational Killing vector (/>", has a family of helical Killing vectors + ^4'°', for each O. Our primary concern, of 
course, is with binary systems, spacetimes for which and 0" are not themselves Killing vectors, although, for one 
value of n, k" =t'^ + ncj)" is. 

We have emphasized that spacetimes with a helical Killing vector cannot be asymptotically flat in the exact theory, 



and a theorem by Gibbons and Stewart |18|, showing that X (null infinity) cannot be periodic, makes this claim 
precise: No spacetime can have a X (and hence no spacetime can be asymptotically flat) if it is vacuum outside a 
compact region and has a helical Killing vector. We can, however, use the Killing vector fc" to define as follows the 
future and past horizon and the future and past domains of outer communication of a spacetime with a helical Killing 
vector. 

Definition II. 2 A point x £ M is in the future (past) domain of outer communication, T>^ if some future- directed 
(past- directed) timelike curve c(A) through P eventually exits and remains outside the history T of each sphere S : 
That is, for each history T that encloses P, there is some Ao for which c(A) is outside of T , all A > Ao. 

Definition II. 3 The future (past) event horizon Ti.^ is the boundary of the future (past) domain of outer communi- 
cation. 



■^Without the requirement on the timelike character of any Killing vector can be written in the form f° + Q(j>" . To restrict 
helical to the vector fields in which we are interested, we had to exclude the spiral Killing vectors of Minkowski space that have 
the form s" + <f)", with s" a constant spacelike vector. 
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Proposition II. 2 Let the history T of a spacelike sphere lie in . Then = dI^{T). 

Proof. Denote by int(T), the points inside a history T . It suffices to show that /~ {T>^) f] int(T) = /~ (T+) p| int(T). 
For any P G /"(!?"'") p| int(T), there is a timehke curve from P that exits T and hence intersects T. Thus P £ 
I-{T +) n int(r); and, from T <Z V+ ^ {T) C I^{T>), the result follows. □ 
The main result of this section is that H.^ are Killing horizons and hence that they satisfy the zeroth law of black- 
hole thermodynamics: That is, the Killing vector /c" is tangent on H.^ to the null generators; and the associated 
surface gravity k, defined by 

k'^Vpk°' = Kk", (2) 

is constant on each connected component of Ti.^. 



To prove that Ti,^ is a Killing horizon (Prop. [1.5 below), we will use an an alog ous theorem proved by Isenberg 
and Moncrief |19 2^ and a strengthened version by Friedrich, Racz, and Wal d (henceforth FR&W), for a class 
of spacetimes with a compact null surface (see also earlier work by Hawking p2[). Following FR&W, we first show 
that the spacetime N, gap covers such a compact spacetime. Although our spacetime is not in the class they study, 
FR&W note that their asymptotic conditions can be relaxed, and we easily extend their proof to spacetimes of the 
kind considered here. 

For convenience in matching our definition and proof to that of FR&W, we consider a subspacetime N = 
/■'"(Af ) Pi intT, for some T that encloses the fluid and black holes. By choosing a future set, we keep all black 
holes but discard the bifurcation horizon and white holes that are part of the full spacetime. (To obtain the corre- 
sponding results for white holes - for the past horizon ~ one exchanges future and past). When the surface gravity 
K is nonzero, the past-directed null generators reach the bifurcation horizon of M in finite affine parameter length. 
This means that in N, they are past geodesically incomplete, and that past incompleteness is one of the conditions 
required for the FR&W proof. The Isenberg-Moncrief version does not require past incompleteness, but does demand 
that the horizon be analytic. N,gaf3 satisfies the following conditions that define a spacetime of type A'. 

Definition II. 4 A smooth spacetime N., g^fj will be said to be of class A' if it has the following properties. 

(i) The spacetime has a Killing vector field k" that is transverse to a Cauchy surface^ 

(ii) N = I+{N) 

(Hi) There is a history T for which N ~ I^{T) 

(iv) The horizon Ti := dI^{N) consists of smooth disconnected components each of which has topology M x . 

(v) The generators of H. are past incomplete. (Alternatively, Ti is analytic.) 

Proposition II. 3 Let N,gap be a spacetime of type A' , satisfying the null energy condition R°'^l°'l^ > 0, all null l" . 
Then on each component of the horizon, there exists a to ^ such that Xto niaps each null geodesic generator of Ti. 
to itself. 

We first need to establish for spacetimes of class A' an analog of Prop. 9.3.1 of Hawking and Ellis showing 
that the shear and divergence of the horizon generators vanish. This implies that the generators are Killing vectors 
of the horizon. Lie-deriving its degenerate 3-metric. 

Lemma II. 1 Let N^gap be a spacetime of class A'. On each component of the horizon, the shear and expansion of 
the null generators vanishes. 



Proof of Lemma Let S* be a Cauchy surface transverse to fc", St = XtiS), and let Bt — Stf]Ti.. Because Xt 

is an isometry, it maps Ti to itself. Then fc" is tangent to Ti, and the family of slices !Ft foliates H. Because Tt is 
mapped to JFc by the isometry Xt'-t, the area of is independent of t. This implies that the divergence 9 of the 
horizon's generators vanishes and that the generators have no past endpoints. Finally, using 9 = and the null energy 
condition, the Raychaudhuri equation (optical scalar equation), 

^=-Raf,l-lP-2aa(3a'^''^^9^ (3) 
implies CTq^ = 0. □ 



Proof of Proposition n.3. Once Lemma [LI is proved, the proof of this proposition is exactly the proof of Prop. 2.1 



in FR&W. □ 



vector field k" is transverse to a hypersurface S if fc" is nowhere zero on 5* and nowhere tangent to S. 
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Definition II. 5 A spacetime N,ga[3 is of class B if it contains a compact oricntahle, smooth null hypersurface Af 
that is generated by closed null geodesies. 



(These causality- violating spacetimes are introduced only as part of the proof of Prop. [1.4; the spacetimes considered 
in this paper as models of physical systems are globally hyperbolic.) 

Proposition II. 4 Let N,gaf3 be a spacetime of type A'. Then (int{N), gap) is a covering spacetime of a spacetime of 
type B. 

The proof is imm ediate: 

Proof of Prop. II. 4 . Because fc" is transverse to a Cauchy surface, Xt has no fixed points for t ^ 0; in particular, 
for t = to of Prop7"^L^, xto has no fixed points. Then the factor space N = mt(N)/xtoi with induced metric gaf3 has 
covering spacetime iiit{N), gap. Because Xto maps each generator of H to itself, H — H/xto is a null hypersurface 
generated by closed null geodesies. □ 



Proposition II. 5 In a spacetime of class A' Ti is a Killing horizon of k". In particular, if, up to a constant scaling, 
/c" is the only Killing vector in N (or in any sub spacetime), then k" is parallel to the null generators ofJi. 

Proof Any neighborhood of a component of the horizon of {\\A{N),gap) that is disjoint from the fluid covers a 
vacuum spacetime of type B. Theorem 4.1 of FR&W implies that in a one-sided neighborhood of that component of 
the horizon, there is a Killing vector K" normal to the the horizon. The puUback if" of to the covering space is 
then a Killing vector on a one-sided neighborhood of the corresponding component of 7i, normal on 7i to Ti: i.e., Ti 
is a Killing horizon. If each neighborhood has fc" as its only Killing vector (up to an overall scale) ,then /c" cx X" on 
each component of 7i, implying that 7i is a Killing horizon with Killing vector /c". □ 

Corollary (0th Law). The surface gravity ni is constant on the ith component of Ti. 

Proof The proof of the zeroth law of event horizons given in Bardeen et al. establishes the result for any Killing 
horizon in a spacetime satisfying the null energy condition. 

The first law is the content of the next section. The second law, that the area of a black hole cannot decrease, 
has meaning here only if one extend the definition of event horizon in a way that requires neither a Killing vector 
nor asymptotic flatness. Black-hole thermodynamics of general spacetimes that are not asymptotically flat has been 
examined previously pq-EM] but the results here appear to be new. n 



III. FIRST LAW FOR SYSTEMS WITH A SINGLE KILLING VECTOR 



We consider spacetimes with black holes and perfect-fluid sources, which have a helical Killing vector or, more gen- 
erally, a single Killing vector that is transverse to a Cauchy surface and timelike on the support of the fluid. Although 
such spacetimes will not, in general, be asymptotically flat, one can obtain a generalized first law of thermodynam- 
ics in terms of a Noether charge Q associated with the Killing vector field and with an action for the perfect-fluid 
spacetime. For spacetimes that are asymptotically flat, the overall scaling of a timelike Killing vector is chosen by 
requiring it to have unit norm at spatial infinity. Here, without asymptotic flatness, the overall scaling cannot be so 
determined. Instead, in our discussion of the first law, the choice of a family of spacetimes will include the choice of 
a Killing vector; but readers should keep in mind that nothing in this section restricts the freedom to choose another 
scaling of the Killing vector for each member of the family of spacetimes. 



"^In particular, in the isolated horizon framework, for a horizon with a single Killing vector, one shows the existence of a 
charge E defined on an isolated horizon for which SE = nSA; in our case this is satisfied by the charge 5Qi = SQli + SQki 
defined on the zth disconnected component of the horizon by Eqs. ( ^9| ) and (^^. Our first law, in contrast, relates this change 
in the black-hole charges to the changes in the Noether charge of a sphere surrounding all black holes and all matter and to 
the changes in the entropy, baryon number, and circulation of the fluid. The existence of such a first law depends precisely on 
what is not assumed in the isolated horizon framework: a globally defined Killing vector. 
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We describe a perfect fluid by its four-velocity and stress tensor 

T''^ ^ eu'^u^ + pq'"^ , (4) 
where p is the fluid's pressure, e its energy density, and 

is the projection orthogonal to u". We assume that the fluid satisfies an equation of state of the form 

P^p{p,s), e^e{p,s), (6) 

with p the baryon-mass density and s the entropy per unit baryon mass. (That is, p :— mBn, with n the number 
density of baryons and niB the average baryon mass.) 
Given a family of perfect-fluid spacetimes specified by 

Q(A) :=[g„0(A),ii"(A),p(A),s(A)], (7) 

one defines the eulerian change in each quantity by 6Q := ^Q(A). 

We introduce a lagrangian displacement ^" in the following way: Let Q Q(A), and let 'I'a be a diffeo mapping 
each trajectory (worldline) of the initial fluid to a corresponding trajectory of the configuration Q(A). Then the 
tangent ^"(-P) to the path A — > 'i>x{P) can be regarded as a vector joining the fluid element at P in the configuration 
Q(A) to a fluid element in a nearby configuration. The lagrangian change in a quantity at A = 0,^ is then given by 

AQ:=^*_aQ(A)U=o = ('5 + £c)Q. (8) 
The first law will be written in terms of integrals over a spacelike hypersurface S, transverse to fc", whose boundary 

^^^SU^B^, (9) 

is the union of black hole boundaries Bi {Bi, is the ith disconnected component of E n 7i^), and a 2-sphere S that 
encloses the fluid and all black holes. Define a scalar t by setting t = on S and requiring k°'V at — 1. 
We can write u" in the form, 

u« = u*(fc" -^w"), (10) 

with u* :— u°'V at and a vector field on S, 

v°'Vat = Q. (11) 
The fact that maps fluid trajectories to fluid trajectories and the normalization u°'Ua = — 1 huply 



Aw" = ^u^w^w^Ag^^. (12) 

One obtains an action for a perfect-fluid spacetime by considering perturbations for which the entropy and baryon 
mass of each fluid element are conserved; and we use this action to define a Noether charge Q associated with fc" , for 
each spacetime Q(A). Then for general perturbations, in which the entropy and baryon mass of each fiuid element are 
unconstrained, we use the charge Q to write a form of the first law for perfect-fluid spacetimes that have one Killing 
vector and a Killing horizon (and that are not, in general, asymptotically flat). 

When the entropy and baryon mass of each fluid element are conserved along the family Q(A), we have 

As = and A(pu"V^) = 0, (13) 

implying 



^The lagrangian change is analogously defined at any Ao: The diffeo = ^'a+Aq^'aq^ maps each fluid trajectory in the 



configuration Q(Ao) to the corresponding trajectory of Q(A + Ao), whence AQ(Ao) := -tt^ -\Q{^ + Ao 



A=o. 
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p ^ 

and the local first law of thermodynamics for the fluid, 

Ae = pTAs + /lAp, 

with 

e + p 



h = 



P 



yields 



Ae Ap 1 . 



e + p p 2 
From these relations, it follows that the scalar density 

1 



is a lagrangian density for a perfect fluid space time. That is, 

SC = ^6iRV^) - A(eV^) + V^ieOV^, 
and (when As = and A{pu°'y/—g) =0), we have. 



^AieV^) = -l{e+p)q'^^Ag^(3 + ^-eg'^^Ag^^ 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



= --T^f'Ag^f^ 

= -^T^PSg^p+^^WpT"^ - V„(T"%). 
That £ is a lagrangian density is then expressed by the equation js^ 

with 



(20) 



(21) 



6" = (e +p)g"^e/3 + TTria'^^g^' - 9''^9''')^p5g^s- 
Ibn 

Now one can associate with C a family of Noethcr charges p^,p^|-p^,p^,p^ H 



(22) 



(23) 



where 



Svr 



(24) 



and i?"(A) is any family of vector fields that satisfies 



^Our Noether formalism is similar to Iyer's extension of the Wald-Iyer work to perfect fluid spacetimes. Like Schutz 

and Sorkin, [||jl4|, however, we use vectors instead of forms, and our lagrangian displacement arises from a map \ from the 
manifold to itself, not, as in Iyer, from a projection onto the manifold of fluid trajectories. 
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-g dX 



(25) 



By choosing i3"(0) = 0, we make Q{X) finite; and, as we will see, Q is independent of the sphere S, as long as 5" 
encloses the fluid and any black holes. Outside the matter. 



(26) 



where, is the covariant derivative of the metric <?q/3(0). 

The generalized first law will be found by evaluating the change SQ in this Noether charge, allowing perturbations 
that change the baryon number and entropy of each fluid element. We restrict the gauge in two ways: We use 
the diffeomorphism gauge freedom to set dk" = 0. The description of fluid perturbations in terms of a lagrangian 
displacement ^" has a second kind of gauge freedom: a class of trivial displacements, including all displacements of 
the form fu", yield no eulerian change in the fluid variables. We use this freedom to set At = 0. Because 6t = {t is 
not dynamical), this is equivalent to the condition ^* — 0. Eq. ( |T^ ) now implies 



—7- = -^"u'^Agais- 



(27) 



Then, from Eqs. (|T|) and (|^, we have Au" = AM*(/fc" + w"), while, by Eq. ([o|), Au" = A[u*(fc° + w°)]; thus 

A(fc"+t;") = 0. (28) 
For perturbations that include changes in baryon number and entropy, Eqs. ( p^ and ( |20| ) are replaced by 



and 



-5 



- ^o"? Ac/q^ H , 



-J—A{pu'^) - l-T'^^Sgo^p + e„V^T"^ - V„(T«%); 
u y—g 2 



and the change in the lagrangian density becomes. 



1 



^ -.5C = -pTAs ~ ^!^A{pu' V^) - -^{G^^ - ^^T'^f)5g^p ~ ^V^T"^ + V^G". 
To find the change 5Q in the Noether charge, we first compute the difference, 



(29) 



(30) 



(31) 



(32) 



between the charge on the sphere S and the sum of the charges on the black holes Bi. As we show below, this quantity 
is invariant under gauge transformations that respect the Killing symmetry. Write Q — Qk + Ql {Qk the Komar 
charge, Ql an additional contribution involving the lagrangian density), with 



J^TT Js Js 



(33) 



(34) 



From the identity 



(35) 



we have 
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fSTT Jq^ 87r 



Now 



= e fe^dS'a + (e + p)u°'uf3V^dSa, 



whence 



(38) 



(39) 



and 



'^(Q/f- VQif,) = - / SCd^x+ f A [(e + p)M"u0z;^d5„] 

-5 ( ^{G^fi-%^T^fi)kPdS^. (40) 
The second term on the right of Eq. ( p^ ) is given by 

= / v^(fc"e'3 - fc'3e")d5„ 

= / VpQ^k'^dSo. - I £ke"dSc., (41) 

where we have used the relation Vak" = to obtain the last equality, and V^jG'^ is given by Eq. (|3l|). Then, adding 
Eqs. (^ ) and (|4l|), and using the relations 

A[{e+p)u°'uf3v'^dSa] = huf3V^A{pu°'dSa) + A{hu0)pu"dSo, + {e + pju'^up £k£.^ dS^, 

where Av^ = —Ak^ = £k£,^ is used and 

£fce"d5„ = (e + p)q"/5 £^^'5^5^ ^ {e+p)u"up£kfdSa,, 

where £k£,^ V/ji = is used, we obtain an expression for S{Q ~J2Qi)- 

i 

S{Q - "3') / pu^dSa + hupv^^^ Aipu^dSa) + v'^A{hup)pu°'dSa 



(42) 



(43) 



-^S I {G''fi~8TTT''p)k^dSo.+ I 



Idtt ^ 



(44) 



We next evaluate the black-hole charges Qi. Recall that, by Prop. (II. 5), fc" is tangent, on each disconnected 
component Tii, to the null generators of the horizon, with surface gravity Ki given by 



(45) 



On each Bi, let n" be the unique null vector field orthogonal to Bi and satisfying ricfc" = —1. The area element of Bi 
is then 



dSai3 = ^{kanp - k/3na)dA. 



(46) 
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Using the Killing equation, V^fc'^ = V'^fc'^l, and Eq. (0) to evaluate the integrand of Qki, we have 



implying 



Finally, following Bardeen et al. ||2j], we show that 

Stt 

Using (S(Vafc^) = ^^{akfi\) = V^afJfc^j], we have 

6k, = 5(n"/c''V„fc/3) 



(47) 



(48) 



(49) 



(50) 



Because the horizon is unchanged in our gauge, and fc^ is parallel to the null normal to TY^, bka — aka, some function 
a on Hi- Then 

6n°'k^V akfj — —Sn^Kika — Kitx^Ska — —Kia 
^ -an°'kfjVak'^ = -xeSkfjVak^ 



(51) 



where, in the last equality, we have used £kSka = 0. From Eqs. (pOj ) and (|5^) and from the vanishing of Saaf3 and 
66, we have 



Now 



SQu = / {k"e^ - kPe'')dSaf3 

^hjs '""^^^'^''^ ~ 9^''9'')^^^9Se \{kan0 - kpn^)dA 

STT 

The first law now follows from Eq. (||) for d{Q - ^ Q,), Eq. ^) for Q^,, and Eq. (||) for QLi 



— KiSAi 
Stt 



(52) 



(53) 



P' f 



— Asu^dSc. + [— + hufsv^ ] A{pu"dSa) + yf^ A{hu(3)pu°'dSa 



Ibn 



(54) 



When the family of spacetimes satisfies the field equations, the last line vanishes and we obtain a first law of thermo- 
dynamics in the form 



/ p^Asu°'dSa + f 4 + huBv'^] A(pu'^dSa) + v'^A(huB)pu°'dSa + V —K^SAi. 



(55) 
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Equivalently, writing 



T:=—, fJ--^— = 1 — ' 56) 



and 

rfMe := pu^dSa, dS := s dMe, dCa ■= huadMs, (57) 

we have 

SQ= I \TAdS + flAdMB + v^AdCJ + V —k,SA^. (58) 
is i 

The relation betwe en thi s form and that for an asymptotically flat spacetime with two KiUing vectors, and 0", 



wiU be found in Sect. IV B 



We noted above that the difference S{Q — Qi) is gauge invariant. In fact, we can see as follows that 6{Qk — J2i Qki) 
and 5{Ql — X^i Qm) are separately invariant under gauge transformations that respect the symmetry fc". The gauge 
transformation associated with a vector field 77" is given by 

5,2 -£„Q, r(^) = -^"- (59) 

The corresponding lagrangian change in any quantity is then identically zero: 

A, = (5^ + = 0. (60) 

From Eq. (pl|), (|3l|), and ( |60| ) the change in 5{Ql — Qu) due to a gauge transformation is given by 

5{Ql - V Ql^) = / V^e^fc^dS^ = / 5^C d^x, (61) 
j is 

when the field equations are satisfied. Decomposing rj in the manner 

Tj" ^Tjf'Wptk'' + fi", (62) 
with fi^S/at — 0, and using £^77" = 0, we have (5,,£ = £,,£ = S/a{Cfi"-), 

5{Ql - V Ql») = / do,{Cr)d'x = 0, (63) 

because £ vanishes outside the fluid (on dYj). 
Similarly, from Eq. (p7|), 

5(Qif - V Qa'O = -^S f R^fik^dS^. (64) 

Again, for a gauge transformation that respect the Killing symmetry, the right side is an integral over the boundary 
dY, of a quantity that vanishes outside the fluid. 

Lastly, we verify the asertion made previously, that Q is independent of the 2-surface S on which it is evaluated, 
if S encloses the fluid and any black holes. This is immediate for Qx from Eq. js^) and (^8|). For Q (and Ql)^ it 
follows from the fact that Q = Qk at A = 0, together with the implication of Eq. (jsj) that ^ = 5Q is independent 
of S along any sequence of equilibria Q(A). 

A. First law in Hamiltonian framework 

In applying the first law to spacetimes that are spatially conformally flat, we will need to write it in a 3+1 form, 
with metric ^ab on S and its conjugate momentum tt'^^ as independent variables. Until Eq. (fz^ ) of this section, the 
vector field fc" that generates time evolution is not assumed to be a Killing vector. 
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Let S = So be a Cauchy surface transverse to k", and let Sj = XtC^), with Xt the family of diffeos generated by 
k". Denote by 7ab(t) the spatial metric on Ej. Let n" be the future-pointing unit normal to this foliation, and recall 
that one can identify spatial tensors on St with spacetime tensors that are orthogonal on all of their indices to n". 
In particular, the projection orthogonal to tIq, 

7a/3 = 5q/3 + '^a"/3: (65) 

is the 4-tensor associated with the family of 3- metrics jab{t) on the slices Sj. Although fc" is not everywhere timelike, 
the fact that it is transverse to a family of spacelike hypersurfaces means that we can introduce a nonvanishing lapse 
a and a shift oj" that relate dt = fc" to in the usual way, 

/fc" = an" + cj", LU°'na = 0. (66) 

Then, in a chart {t,x^} for which St is a i ^constant surface, the metric g^p — 7q/3 — naUp has the form 

ds^ = -a^dt^ + 7y [dx' + oj'dt){dx^ + Jdt). (67) 

With Da the covariant derivative of the spatial metric 7ab, the extrinsic curvature of Sf is given by 

Kab = -■^£n7ab = T7-(-7ah + DaUJb + DbUJa), (68) 

2 2a 

where 7a;, is the puUback to S of £a:7q/3, vanishing when fc" is a Killing vector. 
By taking as independent variables the quantities tt"'', 7a6, a,, with 

^ab ^ _l^j^ab _ ^ab (gg) 

we now generalize the derivation of the first law to permit independent variations of 5tt°-^ , Sjab, 5a, Su!"". 
In terms of Hamiltonian metric variables, the gravitational lagrangian density takes the form ||33| 

= 7r''^7afc - aHo - c^aCg + Da{-2D''aj^ ~ 2w''^% + wV) - n, (70) 

where 

Hg -2G"Vr^;3 7^ = - 'i?7^ + [nabT:'''' - (71) 

eg := -2G"f'^^-nf, 73 = -2DbTT^\ (72) 

Regarding C = {——R — e)\/~5 a function of tt"'', 7cib, a, ixi" and the fluid variables, we rewrite Eq. ( ^ ) in the 
167r 

manner 

(5£ = -a^ipTAs - -^A(pu* a75) 
u 

+ -^{-San - SLu'^Ca + (57r"^[7a6 - DaUJb - DbUOa - 2a{-Kab - ■^7a67r)7"^] 
iDTT 2 

- ^7a6(G''' - 8^5''^)a7n - pT'^^ a^i"- + i?ae"73 - ^(^^-&-^„,)' . (73) 

iOTT 

Here, denoting the puUback to S of a a by ^"'^a, we have set 

PH ■■^T^pn'^nf', -r„^7^"r/, 5,, := T„/37:7f , (74) 

n:=HG + lGTrpH-fK := - IGnfj^ ; (75) 

and the remaining quantities in the last line of ( [73[ ) are given in terms of (tt"'', 7ab, a, w") by 

9" = -^{[ - 2<5(i:>'^a7i/2) + {iv^jbcSn'"' + -kSuj" ~ 2TT\Su;'')h-i 
Ibn 

+ a(e+p)g^e'-a^"/6}, (76) 
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and 



For a Killing vector and A"" any vector field Lie derived by fc", we have the identities 



(77) 

(78) 
(79) 



where 



as Joe 



(80) 



(81) 



and dSa is along the outward normal to dT, in E. In particular, the vector 9a of Eq. ( |76[ ) is related to 0" by 
©a = a0a7a + WaOa^^", implying 



- / Q^dSa. 5QL^ 

Js 

Qk can be expressed in terms of (tt^^, 7a&, a, lj") by writing 



Q'^dSa. 



(82) 



(83) 



with Kab = -(TTab - 57ab7r)7 ^ . Then 



with (34 1 



87r 7s 



R"pkf^na,\^=DaiD''a-K\Lo' 



(84) 



(85) 



We can verify directly that R" pk^Ua takes the form (8E ), when written in Hamiltonian variables, using the Hamil- 
tonian forms already given for Tic, Coa, and G"'^ Eq. M) implies 



Eq. (^) gives 



labG' 



ah 



la. 1 Al a a. 



(86) 



(87) 



and substituting this and the forms ( |7l| ) and (|72|), of Tic and Cca hi Eq. (^6|), we obtain 
i?"/3fc''«a|s = ^a(^"a + ^'•fcC^S-s _ i7rL^''7-3) = Z?„(i?''a - K^J"). 

Consequently, Eq. (|o|) holds with R and G"" pria given by Eqs. (|7^), ([7l|), and ([t^, and with A^ct^ defined as a 
function of (706, a, w"), independent of 7r°^. 

Finally, combining Eq. (^0|), Eq. (^Sj) and Eq. (^) in the lagrangian derivation, as we obtain Eq. ([H^), 
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S^Q-Y^Q^j = J {TAdS + flAdMB + v^AdCa) 

16tt 2 

+ alS-iabiG"^ - Stt^"'') + 167r^„V0r"'3]^5 ]Sx. (89) 

Here the last two integrals in Eq. (||) are combined by using 25[{G°'p - d,TTT°'p)k'^na] = <5aH + aSTi. + (5w"Ca + u^SCa- 
When the field equations are satisfied, and tt*^^ is given by 

tt'^^ = --{D^'^u^^ - -i°'^D^uj'')-i^l'^, (90) 
a 

we have 

SQ^ J [TAdS + p,AdMB + v^AdCa] + J2 '^^^^^ (91) 



IV. APPLICATION TO THE INSPIRALING BINARY BLACK HOLE - NEUTRON STAR SYSTEM 



A. Comparing configurations in quasi-stationary systems 



Our study of a generalized first law was spurred by the fact that equilibria stationary in a rotating frame - spacetimes 
with helical Killing vectors - are used in several approaches to binary inspiral. In each of these cases, one approximates 
the inspiral phase of binary coalescence by an evolutionary path through a sequence of equilibria. The first law has 
a strikingly simple form when used to compare such dynamically related spacetimes: For isentropic fluids, dynamical 
evolution conserves the baryon mass, entropy, and vorticity of each fluid element, and we show that the first law 
becomes 

5Q^^Yk,5A,- (92) 

i 

or 



5Q = 0, (93) 

for perfect fluid spacetime with no black holes. In the gauge that we have chosen ((5fc" = 0), when the spacetime is 
asymptotically flat and fc" has the asymptotic form P + Q.cj)"' , with and timelike and rotational Killing vectors 
of a flat asymptotic metric, we find 

5Q = 5M - Vl5J, (94) 

with M and J the mass and angular momentum at spatial infinity. In particular, the first law in this form describes 
(i) comoving binaries, flows with = 0; and (ii) irrotational binaries, potential flows hu^ — V(j$, with 

A(/lMa) = Vc.A$. (95) 

For an isentropic fluid, conservation of rest mass, entropy, and vorticity have the form 

£m(p\/^)=0, £„s = 0, £„Wa/3 = 0, (96) 

with the relativistic vorticity oJq,^ given by 

= qa^qp^ [V^ {hus) - V5 Qiu^)] = Vq {hup) - (hua) ■ (97) 

The perturbed conservation laws have the first integrals 

A(pm"v^) = 0, As = AtJa/3 = 0, (98) 
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appropriate to the difference between two flows tliat are related by a dynamical evolution. It immediately follows 
that the first and second terms of Eq.(^5|) vanish for isentropic flows. 

To see that the third term vanishes when the perturbed vorticity vanishes, we use [d, = to write 

= AcJap = VaA (hup) - VfjA (hUa) , (99) 



implying Ahua — VaA^, as in Eq. (95). The third term in Eq. ( |55| ) can then be written 

/ v'^A{hup)pu°'dSa = I f^V^A$pu*/fc"(iS'a (100) 

= / Vp{v'^A^pU*)k°'dSo, - / V/3(v^pM*)A$fc"d5a (101) 

The first term in this last equality vanishes, because it is the integral of a total divergence. (Write (V^A^)fc" — 
Wfj{A^k" — A"k^) and use Stokes' theorem; or, more concretely, write k°'dSa = \/~-gd^x ) For the second term, 
recalling the definition of in Eq. (Eoh, we have 



Vp{vPpu') = Vp{pu^) - £k{pu') - pu'Vpk^, (102) 

with each term on the right separately vanishing. 

Thus, for spacetimes related by a perturbation that locally conserves baryon mass, entropy and vorticity, the first 
law has the form (|9^), as claimed. 

B. Asymptotically flat systems 



We will use the 3+1 formalism of Sec. (Ill A) to evaluate 5Q — 5Qk + 5Ql- In the post- Newtonian and in the 
Isenberg- Wilson-Mathews spacetimes that have been used to describe binary systems, the 3-metric has the asymptotic 
form 

7a6 = .fah + 0(r-l), (^Q3) 

where r = (Sijx'^x^Y^'^ , with {x*} a chart for which fij — Sij. 



By writing fc" = an°' + w", as in Sec. (Ill A), we choose a shift u'^ associated with a comoving chart at spatial 
infinity. That is, 

Lu" = rKj)" + (3", where = 0(r"^), (104) 
and 0" is a rotational Killing vector of the flat metric fat 

(j)'' ^x\d2T -x^^diY. (105) 

The extrinsic curvature and lapse have asymptotic behavior 

ifafc = 0(r-'), a = l + 0(r-i), Daa^O{r-^). (106) 

To evaluate 5Q, we first define two asymptotic masses and the asymptotic angular momentum. A mass Mk seen 
by a test particle in Keplerian orbit is associated with the asymptotic form of the lapse, 

Mk — / D^a dSa = lim — / r^dfl, (107) 

where 

/ :— lim 

'Sr 

with 5*^ a sphere of constant r. In terms of the metric potentials, Mk has the form of the Komar mass associated 
with an timelike asymptotic Killing vector t". 
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The ADM mass is computed from the 3-metric: 



= - Urn -!- / drij r'^dn. (108) 



2tt 



FinaUy, the angular momentum associated with the asymptotic rotational Killing vector is given by 

J = -^ f Tr\ 4>''dSa = ^ f K\ cjJ'dSa. (109) 



As in the first equality of Eq. (jsj) , we have 



Qk = / {-D'^a + K''bLo'')dSa = ^ / {D'^a - K\nct>'')dSa, (110) 



whence, by Eqs. (107) and (hO 



Qk = \mk - nj. (Ill) 
We obtain SQl from Eqs. ( ^2| ) and ([76|). Using (jy^-DaV — and the asymptotic behavior given above, we have 

+ (/""/'' - r'')Db5-icd} + 0(r-3), (112) 
,5Qi = -L[_2(5 / D^adSa + ^sj K" b^c/y^'dS a - 2n 5 f K''b(t)''dSa 

IDTT J oo OQ oo 

+ / (/""/"' -r'r')^bW^a] 

'J oo 

= -^5Af/< + s{nj) - nsj + 5M ADM 

= SMadm - ^SMk + Snj. (113) 



Adding Eq. (113) to ^(111), we have 



SQ = ^SMk - S{nj) + SMadm - ^SMk + SflJ 

= SMadm - ^SJ, (114) 

in agreement with the usual first law. 

C. The first law for spacetimes with a conformally fiat spatial geometry 



As mentioned earlier, several groups have recently obtained quasi-equilibrium sequences ]6[-|9[|, approximating bi- 
nary inspiral by a sequence of Isenberg-Wilson-Mathews spacetimes (henceforth IWM spacetimes), spacetimes whose 
3-geometry is conformally flat and whose five metric potentials satisfy a truncated set of five field equations. More 
precisely, the metric of a IWM spacetime satisfies the constraint equations and the spatial trace of the Einstein equa- 
tion, together with the maximal slicing condition for its conformally flat slices; and its matter satisfies the equation 
of motion, V/jT"^ = (see e.g., Isenberg or Flanagan |]36|| ). 

As Detweiler has pointed out, when the spacetime has a helical (or timelike) Killing vector, one cannot in general 
solve all of these equations simultaneously for a metric with conformally flat spacelike slices. One must omit one 
relation to accommodate the new constraint that the existence of a Killing vector imposes on the extrinsic curvature, 
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Kab- We note first that, if one omits the K = condition, the resulting spatially conformally flat spacetime satisfies 
an exact first law, despite the fact that only a truncated set of field equations are imposed. 

In the second part of this section, w e no te that one can alternatively retain the K — condition if one simply 
defines a tensor Kab by the form (Eq.( 123) below) that the extrinsic curvature would take in a spacetime with a 
helical Killing vector foliated hy K — slices. We show that the first law is exact in this framework. This is surprising, 
in view of the artificiality of the definition of Kab and the fact that one component of the Einstein equation is not 
satisfied in the IWM framework. 

In each case, one has a spacetime foliated by hypersurfaces whose spatial metrics have the form 

7ab = V'^'/ab, (115) 

with fab a fiat metric. The corresponding 4-tensors, 

lap = -0*^/0/3, (116) 

are Lie derived by the Killing vector fc": 

£/c7a/3 = 0, £fcV = 0, £kfaf3 = 0. (117) 

In particular (although we will not use the fact in this section), 

A:"=t" + f70", (118) 

with ((>"■ a rotational Killing vector of fab- 

In the first case (with K not required to vanish), the spacetime satisfies on each St the equations 

H^O, Ca^O, 7a6(G'''-87rr''^) =0, VpT'''' = 0, (119) 

together with the relation ( pO| ) expressing tt"'' in terms of the metric. Because 

Slab = 4^7a6, (120) 

it is exactly this set of equations that occur in the action and in the 3+1 form of the first law (^), when one compares 
two spatially conformally fiat spacetimes. 

Finally, comparing asymptotically flat spacetimes of this kind, with no local change in entropy, baryon number, or 
circulation, we have 

SM = n6J + J2'^iSAi- (121) 



We consider next solutions [tt"'^, jab, Oi, w", e, u"), to the same set (119) of equations, now with tt = 0: 

^ah^_^-afc^i^ (122) 

with Kab the tracefree part of the extrinsic curvature: 

1 2 

Kab = 7^(i5aCJ6 + DbUJa - -JabDcUj''). (123) 

One writes H, Ca, 7afc(G"'' - SttT"''), and V^jT"^ as they occur in the Hamiltonian formalism (for the metric), as 
fimctions of (tt"'', 7af,, a, a;°) and the matter variables; one substitutes for -fab and 7r°^ the expressions 

lab = i^* fab, 7r'''' = -^{D-u;' + D'>u;--l-f-'D,u;'')-f-i; (124) 
and one solves the resulting system of equations for (-0, a, w"; e, w"). 

n^O, 4 = 0, (G'^' - SnS'^'hab ^ 0, VpT"^ = 0, (125) 
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where 



n = n(TT ,jab,a,i^ j, L =C (tt ,7^6, a, ;e,u j, G 7q6 = G 7afc(7r ,7a6,a,a; ;e,u j. (lib) 

Then, for a family of such solutfons, the quantities 6a, Slu, and S'^ab — 4:^lab occurring on the right of the first law 
( p9| ) multiply expressions that vanish. Because Jtt'^^ is traceless, the expression involving Sn"-'' has the form 



1 



167r 

Eq. (p9|) thus yields 



5TT''''[DaLUb + DbUJa + 2a"f ^TTab]^ --^6TT'"'[DaLUb + DbUJa- -JabDcUj" + 2aj 2 TTafe] = 0. 



16n 



SQl-^S J R°'0kf^dSc = J [TAdS + -pAdMB + v^'AdCa]. 



(127) 



(128) 



To recover the first law in the form 



we must show that 



f ITpkl'dS^ - -^5 f V'^kl'dSo.f}. 



(129) 



(130) 



This is not obvious, because, in replacing the extrinsic curvature by its tracefree part, we invalidate the Killing identity 



(131) 



Remarkably, however, the Ua components of the two sides of this inequality differ by a divergence; and the asymptotic 
behavior of the spacetime implies the equality 



1 

"8^ 



Wkf^dSc^f 



1 

'8^ 



ik^pk^dSa 



(132) 



That is, from Eq. 



and 



Then 



we have 

Wl3V"k>^nc 



As noted in Sect. (IVB), 



ITak'^na 



DaiD''a-K\u;''). 



uj" ^ n(t)'' + (3" , with P^^Oir-^), 
where 0° is a rotational Killing vector of the flat metric fab', and K — 0(r^^). |^ We then have 



(133) 



(134) 



(135) 



(136) 



^ If, however, one allows a nonzero 3-momentum, with boosted-Schwarzschild asymptotics, then ^ = 1 + f{f)/r + 0{r ^) and 

and /3" = 0(r~^). Because K is given by —( — Q(l>'^Datp + DaP"^), one can only demand K = 0(r~^), /J^if = 0(r^^), allowing 

a ip 

a finite contribution to ■ 
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^TT j ' Stt 7s 247r Jq^^ 



1 

as claimed. 



s 



R°'fskPdSa, (137) 



From Eq. (lOS), we conclude 



5M = nSJ, (138) 
along a family of conformally flat solutions to the IWM equations, written in terms of (tt"'', -0, a, w°). 

Note that the equation ^pV aT°'^ = is satisfied, because, for an iscntropic fluid, the equation of hydrostatic 
equilibrium, conservation of rest-mass, and the one-parameter equation of state together imply Vq.T"'' = 0. To see 
this explicitly, we decompose the divergence of the stress tensor as follows: 

V„T"^ = {q,3-y - Ufsu^)VaT°'^ = P [v?V.y{huf3) + Vfjh] + uphV^ipu") - pTV p s (139) 

In constructing an isentropic (s = const) equilibrium model, conservation of rest mass is assumed, and a barotropic 
equation of state p = p{p) is used. Helical symmetry and the assumption that the fluid flow is either co-rotational 
or irrotational then leads to a first integral of the Euler equation u'V ^{huf}) + Vph — 0. It is this first integral, 
specialized to a conformally flat metric, that is solved in the IWM formalism, implying that Vq.T"'' = 0. Thus, as 
claimed, all terms involving the field equations vanish in Eq. (jsj), and the first law holds for IWM spacetimes in the 
form (KF 



As in the exact theory, when the system includes black holes, the KiSAi terms refer to Killing horizons. The IWM 
spacetimes do not satisfy the Raychaudhuri equation for the null generators of the horizon; as a result, as noted in 
the introduction. Killing horizons in IWM spacetimes need not co-rotate with the orbital motion. 



V. DISCUSSION 



The first law can be used to deduce a criterion for orbital stability for the asymptotically fiat models of binary 
equilibria discussed above, using a theorem of Sorkin ^7\ . Consider a one-parameter family Q(A) of binary equilibrium 
models along which baryon number, entropy and circulation are locally constant (the lagrangian changes As, AdMs, 
and AdCa vanish). Suppose that J = at a point Aq along the sequence, and that fiJ 7^ at Aq. Then the part of 
the sequence for which flJ > is unstable for A near Aq. 



The result relies on a first law in the form 

dM = ndJ (140) 

and on the fact that the equilibria are extrema of mass with J constant. As we have seen, this is the case for a config- 
uration space in which baryon number, entropy, and circulation are fixed for each fluid element. For asymptotically 
flat models with one or more black holes, if one also fixes the area of the horizon along a sequence, then the same 
criterion above can be used to diagnose stability. 



In general, the proof of the theorem shows only that the spacetime is secularly unstable on one side of the turning 
point. In the present context, however, the theorem shows the existence of nearby configurations with lower mass 
that can be reached by perturbations that conserve baryon number, entropy, and circulation; this suggests that the 
criteria locates the onset of dynamical instability. 



When one models stationary binary systems in full GR, the lack of asymptotic flatness leads to several ambiguities. 
For binary charges in Minkowski space, one can obtain a one-parameter family of equilibria if one simply replaces 
asymptotic regularity (finite energy) by a condition that the electromagnetic field be given by the half-advanced -I- 
half-retarded Green's function. In GR, it remains to be seen whether one can find an analogous asymptotic condition. 
Simply requiring equal amounts of ingoing and outgoing radiation is a weaker condition even in Minkowski space; in 
GR one must have asymptotic conditions as restrictive as asymptotic flatness to avoid ambiguity in each asymptotic 
multipole. Finally, as mentioned in Sect. Ill, the helical Killing vector has an arbitrary scaling that one must resolve 
to obtain a unique value for the charge Q of the first law. 
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APPENDIX A: VIRIAL RELATION IN IWM SPACETIMES 

In this appendix, we derive a virial relation for quasiequilibrium states in IWM spacetimes. Incidentally, we show 
that the virial relat ion is equivalent to the relation Mk = Madm- 

As described in IV C, we use a 3+1 formalism, with 3-metric 7ab = ip'^fab, and with a helical Killing vector that 
has the form 



(Al) 



where 0" is a rotational Killing vector of the flat metric fat- Throughout this appendix, we use Cartesian coordinates 
t, {x*} for which fab has components fij = 6ij and = 5q. In the IWM formalism, one imposes the maximal slicing 
condition if = on the family of t =const surfaces St; and, instead of solving the full Einstein equation, one solves 
the Hamiltonian constraint, the momentum constraint and the equation for the slicing condition. Here, however, as 
in Sect. [VC, to obtain a set of equations consistent with the existence of a helical Killing vector, we replace the 



extrinsic curvature in this set of equations by its tracefree part, K^^ . 
The basic equations are then 



7 



Ax = 2^xr{pH + 28^^) + -xV"X/if/ = S; 



(A2) 
(A3) 
(A4) 



where A denotes the flat Laplacian for three space, — det(7y ) — tp^, and x = Q^V'- (See Eqs. ( [t^ ) for definition of 
Ph, ji and Sij.) The energy-momentum tensor is assumed to be nonzero only inside the light cylinder {x^ + y^Y^"^ < 



The shift vector of Eq. (|104D satisfies 

p = —an + 1 . 
The Cartesian components Kf are given in terms of /3* by 



The asymptotic behavior of geometric variables is that of Eqs. (103)- (106) 



V; = l + 0(r-i), 
X-l + 0(r-i), 
/3^ = 0(r-2), 

K,' = 0(r-3), 

appropriate for an asymptotically flat spacetime in a chart for which the total ADM 3-momentum vanishes: 



(A5) 



(A6) 



(A7) 
(A8) 
(A9) 
(AlO) 



®The relation Mk = Madm for stationary and asymptotically flat spacetimes has been proven by Beig js^], and by Ashtekar 
and Magnon- Ashtekar A virial relation relying on this has been derived by Gourgoulhon and Bonazzola po| . 
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(All) 



where daj = {Djr)r^dfl and is computed in Cartesian coordinates. 

We now derive the virial relation and show the equivalence of Komar and ADM mass for quasiequilibria of two 
neutron stars. For r oo, x and ip behave as 1 — M^/2r + 0{r^^) and 1 + A/adm/2?' + 0(?'~^). From this asymptotic 
behavior, we can define M-^ and Madm by the surface integrals, 



M- 



Mabm = — 



1 

2^ 



(A12) 



Since — —Madm + 2Mx, our goal is to show — Madm- 

Qncr 

{'ii^S^+5'^d,xdj'>P)<i^x, 



Using Gauss's law, they can be rewritten in the manner 

1 



^ 27r 
1 



Madm = {xS^ - d,i,d,x)d''x. 

We next derive a relation that will be used several times in the calculations that follow. From 
xi^^K/K/ = y^k^^dj(3\ we have 



(A13) 
(A14) 



XrK^'K^'d'x = j ^K^'djp'd'x 
= - J d,{^k^')P'd^x + j ^k^' P'da, 
= -Svr / ^3^13'd^x, 



(A15) 



where we use the asymptotic behaviors at r ^ cxd and Eq. (A3) to obtain the last line. § without specification of 
a surface denotes a surface integral over dYi: § — §^^^- From the vanishing of the total ADM 3-momentum (more 
precisely, from the vanishing of J^^^ kf^dSj) and from the momentum constraint (|A3|), we have 

(A16) 



= y 3i\/i<i^x, 

which may be interpreted as the linear momentum of a neutron star. 

Using Eqs. ( A13| ) and (A14), we write the difference between Madm and M^ in the form 



M^ - Madm = - 

TT 



= 2 



l^xi^'^S^ + -x'^^'K^'K.' + b''d,^d,x 



d^a 



Vlijkv'' + 3aP} - iv7Ji/3^' + ^5''d4>d,x 



d^x, 



(A17) 



where we use S/J' = jkUi"/''^ /{au*) + 3P, — u*'7"*|_|i = u*'{v^ + (3^) and Eq. (A15). In the following we show 

that the relation M^ = Madm is equivalent to the virial relation. 

To derive the virial relation, we first write the general relativistic Euler equation ^'^^SJ ^T^^ = in the form 

dt{jk^) + dk{3i^v') + dk{a^P) + PHi^'dkX - {rh + 2Si')x^^dkil^ - ^IJidkP' + \xi^S,jdkf' = 0, (A18) 



^The above definition for v'^ is used only in this appendix. Note that v'^ was differently used for spatial velocity vector in 
co-moving frame with fc" as defined in Eq.(fLO|) in main sections. 
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where 7*^ = ^^ip'^. Equation (A18) is a fully general relativistic expression. In the IWM spacetimes, 7'^ = 7.^ = Sij 
and consequently, the last term in Eq. (AlS) is neglected. 

In the following calculation, we choose the x^-axis so that, on some time-slice Ej, it lies along the centers of the 
two members of the binary system. 

As in the Newtonian case, the virial relation can be derived by taking inner product with x'^ and by performing an 
integral over three space, i.e.. 



d^xx'' 



0. 



Below, we shall carry out integrals separately. For simplicity, we omit hats ( " ) on indices in the following. 
(1) First term: Since we assume the existence of the helical Killing vector, we have a relation 

dtJk = -^[80 Jk) + JidkA, 



(A19) 



(A20) 



where we use dicj)^ — and Jk = 3k\Jl- In the present coordinates, (fi 
parts, we obtain 



X , x — ^,0). After an integration by 



x^'dtJkd^x = ^ I {cj,''Jk - Jix''dk^')d^x. 



From a relation x^dk4>^ =(/>'+ 6^^h, we immediately find 

j x''dt{jkVi)d''x = -bn J j2Vid^x = -bnp^s, 

where Pns is interpreted as the linear momentum of a neutron star [see Eq. ( A16| )]. 
(2) Second and third terms: An integration by parts immediately yields 



x^dj{jkv^^)d^x = - y jkv'^y^d^x, 
J x''dk{a^P)d\ J a^Pd^x. 



(3) Fourth and fifth terms: Using Eqs. (A2) and (A4), we can rewrite these terms as 



PH^'dkX - {PH + 25/)xV'afcV = 



AipdkX + Axdki> 



167r 



Taking into account an identity, 

J[{x''dki^)Ax+{x''dkx)Ai^]d^x^ J S'^d.xdjH^'x, 

we find 



1 J^s^ 



2tt 



5'^d^xdjip 



IGtt 



d'^x. 



(4) Sixth term: 



^j^x'^dkP'd'x 



dj{^K^')x^du(3'd'x 



^ / [^/lK^'x'^^k^,^3' + ^K^'d,p']d^x ^K^' x^dkP'da, 



1 

'8^ 



[dk{^/lK.^')x''d,(3' + 2^K^'d,f3']d^x + — j> ^{K^'x'' - K^''x^)dj(i'dak 



(A21) 



(A22) 



(A23) 
(A24) 

(A25) 
(A26) 



(A27) 
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1 

8^ 
1 

16^ 



2^/7 



1 

16^ 



+ ^ iaij^'K^'K/x'^dak - ^ I ^K^'^x^ dj/3'dak, 



(A28) 



where we us e Eqs. (A3) a nd (ATs] ). Because of the asymptotic behavior, the surface terms at r ^ oo in the last hne 
of Eq. (A28) vanish. Eq. (A16) imphes that the center of mass of the system does not move in the x'^ direction, that 
the sum of the momenta of the neutron stars vanishes. 
Gathering the results of (1)~(6), we obtain the relation 



: 



(A29) 



This is the vi rial r elation for a neutron star bi nary s ystem in quasiequilibrium. 
From Eq. ( A17 ), the right-hand side of Eq. ( A29 ) is written as 



= - 



— Mabm 



(A30) 



implying Madm = = Mk, if the virial relation holds. 

APPENDIX B: THE FIRST LAW FOR NEWTONIAN BINARY SYSTEMS 

In this appendix, we derive a first law of thermodynamics for Newtonian gravity. We start with a first-order 
perturbation of the energy of a perfect-fluid. 



E = T + W + U, 



where 



V,$nV,;$n ) dV, I pudV. 



and $N and u denote the Newtonian potential and specific internal energy^. An integral equation 

6 f fpdV^ f AfpdV+ f fA{pdV), 
Jv Jv Jv 

is satisfied for a perturbation. 

The perturbation of the kinetic energy T can be expressed as follows 



(Bl) 



(B2) 



(B3) 



ST 



pv'Av,+ [ lv^A{pdV)+ [ [^^vj'^i{pv')+ p^^v'\/,Vj]dV - (f pv'vjS^^dS^ 
V Jy 2 Jy Jqy 



(B4) 



where we used the relation, 



^°We use for the fluid velocity vector in the inertial frame in this appendix. Note that it was differently used for spatial 
velocity vector in co-moving frame with k" as defined in Eq.(p^ in main sections. 
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2 J ^ 



The perturbation of the gravitational potential energy becomes 

5W= f $NA(pdF)+ / pCV.^Ndl^- / (V^^N - 47rGp) + / V,$N'5$Nd5', 

Jv Jv 47rG Jy AttG Jqv 



'V Jv 
The perturbation of the internal energy becomes 



SU= [ pTAsdV+ [ (u+-) A{pdV)+ [ C^^PdV-(f CPdS^. 

Jv Jv \ P J Jv JdV 



where we used a relation 



as well as a local thermodynamic relation, 



Au = TAs + -^Ap. 

p^ 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



Surface integrals appeared in expressions for ST, SW and 6U are all vanish. Combining Eqs.(B4), (B6) and (B7), we 
have a perturbation of the Newtonian energy integral: 



SE ^ ST + 5W + 6U 



/ pTAsdV^ 


Li 


Jv 





1 o P 

- + $N + W + — 

2 p 



^ A{pdV) + J pv'Av.dV 



^'vjV,{pv') + ep ( v'^^jv, + iv,P + V,$N ^ - ^(V'*N - 4^Gp) 5$, 



= j pTAsdV + J (^^v^ + +U+ A{pdV)+ J pv'Av^dV 

1^ + V,(pz;') dV + f p + i;^ V,i;, + ^ V,P + V,$n ) dV 



47rGjy Jy 

Next, we derive a variation of the total angular momentum J defined by 



J- / pv.c^'dV, 
Jv 

where 0* is a generator of rotation with Cartesian components (jf = (— j/,x,0). The variation of J is 



(5J= / pAv,(l)'dV+ / pw,A</)'dF+ / v.cj)' A{pdV). 
Jv Jv Jv 

Using a relation 

A0' = ,50' + = -£0e, 
the second term of Eq.( B12| ) is rewritten as follows: 



(BIO) 



(Bll) 



(B12) 



(B13) 



pviA(t)'dV 



pv,£^CdV^- / £4pv^e)dV+ / e£dpv^)dV 



pv,e4^dS,+ / e^<k{p^^)dV, 
av Jv 



(B14) 



where we used Vj0-' = 0. Discarding the surface term in the above expression and substituting in Eq.(B12), we have 
a variation of the total angular momentum bJ as follows: 
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6J= / v,(j)' A{pdV)+ / p(j)'Av,dV+ / C^4pv^)dV, 



(B15) 



Finally we write down a general expression for the combination of SE and Q6J, where is a constant parameter, 

6E - nsj = 

pTAsdV + J Qu^ + $N + M+ ^ - A(p(iT/) + J p{v' -n4>') Av.dV 



dp 



dpvi 



at 



dV + 
dV, 



dt 



v^VjV, + -ViP + V,$N 1 dV - 



47rG Jy 



(B16) 



As an application of the above general expression, consider a Newtonian binary star system in circular orbit pi] . 
In this case, the fluid variables Q admit a helical symmetry, namely. 



2 = 0, 



(B17) 



that is, the last integral in Eq.( B16| ) vanishes. When a Mass conservation equation, the Euler equation and the 
Poisson equation for the Newtonian gravity are satisfied, namely. 



-^+V,(pt;') = 0, _5Jl+t,JV,«, = — V,P-V,$N, and V^^n = 47rGp, 
ot dt p 



Eq.(B16) takes a simpler form, 



/ pTAsdV^ 




Iv 





^ + $N + 7i + - - t-i f]0M A{pdV) + I p{v' - n 4>') Av, dV. 



.2 P 

If we further assume that the perturbed flow is isentropic and mass conserving, 

As = and A{pdV) = 
and that the vorticity of each fluid elements is conserved, 

Alu.j = A{VjV^ - V^Vj) = V■iAv^ - V^Av-i = 0, 



then Eq.(B19) reduces to 



6E^Q 5 J. 



Here we have used Eq.( B21 ) to introduce a function for which 



(B18) 

(B19) 

(B20) 
(B21) 
(B22) 
(B23) 



this form of Avi, together with helical symmetry imply that the last term in Eq.(Bl£) vanishes: 

j p{v'-n(t>')Av,dV^J p{v' -n<p')\J,-^ dV j p{v' -il(j)')-^dS, + J (^^ + £n4,p^ ■i'dV^O (B24) 
where we used mass conservation equation and V^^* = 0. 
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